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"*^. '. Abstract 

"^ . In 1853 Sylvester stated and proved an elegant formula that expresses the 

■ ' polynomial subresultants in terms of the roots of the input polynomials. 

Sylvester's formula was also recently proved by Lascoux and Pragacz by using 
multi-Schur functions and divided differences. In this paper, we provide an 

^N I elementary proof that uses only basic properties of matrix multiplication and 

QQ ' Vandermonde determinants. 

o 

^ : 1. Introduction 

^ ■ Subresultants play a fundamental role in Computer Algebra and Computational 

c3 ■ Algebraic Geometry (for instance, see (|a|3|;y;|l3lld;y;llJ;ll2;ll|))- In (|l6l ) Sylvester 

p ! stated and proved an elegant formula that expresses the polynomial subresultants 

^ I of two polynomials in terms of their roots, the so-called double-sum formula. This 

identity was proved also by Lascoux and Pragacz in ()l3h . by using the theory of 
multi-Schur functions and divided differences. 

In this paper we provide a new and elementary proof that uses only the basic 
properties of matrix multiplication and Vandermonde determinants. As apparent in 
our proof, Sylvester's double-sum formula is only one simple step further a particular 
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case, the so-called single-sum formula. Such connection between the single and the 
double-sum formulae was originally thought to be unlikely, as remarked in page 691 
of (|l3t). There have been various proofs for the single-sum formula (|l|; |2; \^ |ll|; |8|). 

The matrix multiplication technique, presented in this papers, has proven to be 
quite powerful in that it is easily generalizable to multivariate polynomials: similar 
techniques were successfully applied to obtain expressions for multivariate subresul- 
tants in roots in (|7|), and the generalization of Sylvester's single and double-sum 
formulae to the multivariate case is the subject of ongoing research. 



2. Review of Sylvester's Double Sum for Subresultants 

Let / = ttrnX"^ + ■ ■ ■ + Oo and g = hnX^ -|- ■ • ■ -|- 6o, be two polynomials with coefficients 
in a commutative ring. The d-th subresultant polynomial SYeSd{f,g) is defined for 
< d < min{m, n} or, ii m ^ n holds, for d = min{?77,,?7,}, as the following 
determinant: 



STesdif,g) := det 
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where ai = bi = ior i < 0. 

By developing this determinant by the last column, it is clear that Sresrf(/, (7) is 
a polynomial combination of / and g. It is also a classic fact that Siesdif^g) is a 
polynomial of degree bounded by d, since it coincides with the determinant of the 
matrix obtained by replacing the last column Cm+n-2d by 



L/„ 



m-\-n—2d 



^m+n-2d ~ ^ ^m+n-2d-l 



X 



m+n— d— 1 



Ci. 



Now, let A = (. . . , a, . . .) and B = (...,/?,...) be finite lists (ordered sets) of 
distinct indeterminates. In (jlGj) Sylvester introduced for < p < \A\,0 < q < \B\ 
the following double-sum expression in A and B: 



\A'\=p,\B'\=q 



where 



R{X,Y):= n (x-y), Rix,Y) := l[{x - y). 

xeX,yeY y€Y 
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In f)16[ ) Sylvester gave the following elegant formula that expresses the subresul- 
tants in terms of the double-sum, that is, in terms of the roots of / and g. 



Theorem 1 (Sylvester's double-sum formula): Let f,g be the monic poly- 
nomials 

f=l[{x-a), g=l[{x-P) eZ[aeA,peB][x], 

where \A\ = m and \B\ = n. Let p,q > be such that d := p + q < min{m, ri} or 
d = min{m, n} if m ^ n holds. Then 

Sres,(/,^) = ^-^ SjW^'\A,B-x). 

\pj 

When p = d and q = 0, the above expression immediately simplifies to the single- 
sum formula: 

Sres.(/,^)= E^(^'^')|^^|;||. (2) 

\A'f=d 

Complete proofs of Sylvester's double-sum can be found in (|l6l : Il3l ). while the 
single-sum formula has various proofs, (|l|; 13 la lll|;|8|). Here we present in Section HI 
an alternative elementary proof for both results. 



3. Notations 

We recall that < d < min{?77,, n} or d 
and Mg denote the following matrices: 



min{?7i, n} ii m ^ n holds. We let Mj 
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We now define 



Sd-.-- 
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Finally, we define for a polynomial p{t) and two lists, T := (71, . . . ,7n) of scalars 
and E := (ei, . . . , e^,) of non- negative integers, the (not-necessarily square) matrix 
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of size V X u: 



{p{t),r)E:- 



iTpiii) ■ ■ ■ luPiriu) 

iTp{1i) ■■■ l"uP{lu) 



For instance, under this notation, if we take i? := (0, . . . , m— 1), we have the following 
equality for the Vandermonde determinant V(r) associated to F: 

V(F):=|(7r^)i<,,<„| = |(l,F)^|. 

When E is of the form i? = (0, . . . , v — 1), we directly write {pit), F)^,. 

We mention the following useful equalities that hold since m + n — d > max(m, n): 



Mg ■ (1, T)m+n-d 
Mt-x ■ (l,F)m+n--d 



{f{t),TU-a 

{git), T)ra^d 
{t-X,T)d 



4. The Proof 



The proof is divided into a series of lemmas which are interesting on their own. 
For an easier understanding, we recommend not to pay attention to signs in a first 
approach. 

Lemma 1: Under the previous assumptions and notations, we have 

SreSrf(/,(7) = (-l)'^+("-'^)(™-^)|5rf|. 

Proof: We denote by Cj the i-th column of the matrix Sd and we replace its first 
column Ci by C[ := Ci + XC2 + . . . + x"^~^"'~'^~^Cm+n-d- This operation does not 
change the determinant of this matrix, and 



C[:- 




We now perform a Laplace expansion of the determinant of the new matrix over 
the first d rows, and we observe that only one block survives, which corresponds to 
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columns 2 to rf + 1 of Mt-x- Moreover, this block is lower triangular with diagonal 
entries 1. Thus 



\Sd\ = (-l)^det 
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= (_l)'^+(n-<i)(--d)sresrf(/,g), 

since the matrix in the right-hand side above is the matrix of (P) viewed backward. 

D 

For simplicity, from now on, we assume / and g to be the monic polynomials 

/ = naGA(^ ~ ck)) 9 = rifles (^ " P) where A and B are hsts with |y4| = m and 
\B\ = n. (As pointed out by a referee, under this assumption one has in the language 
of multi-Schur functions: \Sd\ = 5'id.(m_^)n-d.om-d(— x, — A, — S) (see (|l3ri).) 

The lemmas below generalize in an obvious manner to non-monic polynomials. 
The first one corresponds to Th. 3 in (jllh. We prove it here with a different technique 
that follows from Lemma ^ 

Lemma 2: (Hong's suhresultant in roots (li, Th. 3.1)) 
Under the previous notations, we have 



Sres,(/,^)V(A) = det 
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Proof: We note that \Sd\ V{A) is the determinant of the following product of ma- 
trices: 
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Since {f{t),A)n-d= [«}-'/(«.)] = [0]- 
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By permuting the rows of the second block with those of the third, we obtain 
SreSd{f,g)V{A) = (-l)'^+("^-'^)("-^)|5,| V(A) 
= (-l)'^det 



{t-x,A), 



(git), A) 



m~d 



|m;i 



det 



{x-t,A)d 
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since M'r is a lower triangular matrix with diagonal entries a^ = 1. 
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Let us remark here that the Poisson product formula Res(/, (?) = naGA5'('^) ^^ ^ 
direct consequence of the previous Lemma for the case (i = 0. 

For 5* C T finite lists, let sg (S", T) := (—1)'^ where a is the number of transposi- 
tions needed to take T to S' U {T\S). Here, "U" stands for list concatenation and 
"\" means list subtraction. 

Lemma 3: Let P and Q he two disjoint suhlists of E := {0, . . . ,d — 1) that satisfy 
PUQ = E, and let p := \P\, q := \Q\. Then 



Sresdif,g)ViA)ViB) = (-!)«+(— ^)" sg (P,i5;) det 
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Proof: Recalling that V{B) = |(l,i?)„|, we have by Lemma|21 



Siesd{f,g)V{A)ViB) = det 
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since Mg ■ (1, B)m+n-d = {g{i)-,B)m-d = [0]. Now, since the first matrix is lower triangular 
with diagonal entries 1, we have 



Sresd(/,5) V(^) V(B) = (-l)(— '^)« det 



{x-t,A)d 





(l,^)m+n-d 


(X,B)m+n~d 
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Finally, recalling that {x — t,A)(i 

\i) Ajm+n—d - 

rows yield 



a* ^x — a* 



and 



°r' 
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, the obvious subtractions and permutations of 



Sres,(/,<7)V(^)V(S) = (-1 
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The lemma follows by moving (—1)'' out of the determinant. 

We will also need in the proof the following observation: 
Observation 1: LetT ■= (ji, . . . ,'jd)- Then 

\{x-t,T)d\=Rix,T)\{l,T)d\. 
Proof: The claim follows from 



/ x -71 



Td 



■7d \ / 1 



■^i J V 7?- 



1 X / x'-7i 



n 



(5) 



x-id 



D 



4.1. Proof of Theorem ^ 

For any P and Q disjoint sublists oi E := (0, . . . , d — 1) that satisfy P U Q = E, 
with |P| = p and |(5| = q, a Laplace expansion over the first d rows in Identity (jHl) 
gives that Sresd(/, g) V (A) V (B) equals 

i-irsgiPE) Yl sgiA',A)sgiB',B).\{x-t,A')p\.\{x-t,B')Q\.V{A\A'uB\B') 

A'cA,B'cB 

\A'\=p,\B'\=q 

where a := q + {m — d)n + (m — p)q = {m — d){n — q) (mod 2). Adding over all 
such choices oi P C E with \P\ = p, we deduce that Sresd(/, g) V (A) V (B) equals 

-^ Y, (-1)" sg (P, i?) ^ sg(^'. A) sg{B', By\{x-t, A')p\-\{x-t, B')q\-V {A\A' U B\B') 

\p) P A',B' 
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= ^-^Y.^s{A',A)sg{B\B)V{A\A'uB\B')(Y.^g{P,E)\{x-t,A')p\-\{x-t,B')Q\\ 
\p) A',B' \ P ) 

We observe now that, by another Laplace expansion and Identity (jSJ, 

Yj^g{P,E)\{x-t,J^)p\-\l,x-t,B')Q\ = \{x-t,J^\JB')^\=R{x,J^)R{x,B'%\,A\jB')^ 
p 

Recalling that \{l,A'UB')d\=V{A'UB'), this gives 

,, (f . i-^rS^r,t A'^m o/^ sgjA', A) sgjB ', B)V {A\A' U B\B') V [A' U B') 
Sresrf(/, g) = . 2^ R{x, A )R{x, B ) 



_ i-in-ir ^ ^, ^^.p, ^. R{A',B')R{A\^^B\B') 

where r = {m—p){n—q)+pq—{m—p)p—{n — q)q = {m—d){n—d) since for any finite 
lists X, y, one has V(Xuy) = V{X)V{Y)R{Y,X) = {~iy^^-^^^V{X)V{Y)R{X,Y). 
The claim follows now from the fact that (m — d){n — q) + {m — d){n — d) = (m — d)p 
(mod 2). D 

As a final remark, we mention that if in the previous proof we start with a Laplace 
expansion over the first d rows in Identity Q instead of Identity Q, we obtain in 
the same manner Sylvester's single sum formulation Q. 

Acknowledgements. We are grateful to the anonymous referees for their careful 
reading of our preliminary manuscript and their very precise indications to improve 
our presentation. 
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